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ABSTRACT 

Context. To investigate the disk formation and jet launching mechanism in protostars is crucial to comprehend the earliest stages of 
star and planet formation. 

Aims. We aim to constrain the physical and dynamical properties of the molecular jet and the disk of the HH 212 protostellar system 
at unprecedented angular scales exploiting the capahilities of the Atacama Large Millimeter Array (ALMA). 

Methods. ALMA observations of HH 212 in emission lines from sulfur-bearing molecules, SO 9g - 87, SO lOn - lOio, SO2 82,6 -'ii.i, 
are compared with simultaneous CO 3-2, SiO 8-7 data. The molecules column density and abundance are estimated using simple 
radiative transfer models. 

Results. SO 95-87 and SO2 82,6-71,7 show broad velocity profiles. At systemic velocity they probe the circumstellar gas and the cavity 
walls. Going from low to high blue-/red-shifted velocities the emission traces the wide-angle outflow and the fast (~ 100-200 km s“*) 
and collimated (~ 90 AU) molecular jet revealing the inner knots with timescales < 50 years. The jet transports a mass loss rate 
> 0.2 - 2 X 10“* Mq yr“*, implying high ejection efficiency (> 0.03 - 0.3). The SO and SO2 abundances in the jet are ~ 10“’ - 10“^. 
SO lOii - lOio emission is compact and shows small-scale velocity gradients indicating that it originates partly from the rotating disk 
previously seen in HCO^ and C*’0, and partly from the base of the jet. The disk mass is > 0.002 - 0.013 Mq and the SO abundance 
in the disk is ~ 10“* - 10“’. 

Conclusions. SO and SO2 are effective tracers of the molecular jet in the inner few hundreds AU from the protostar. Their abundances 
indicate that 1% - 40% of sulfur is in SO and SO2 due to shocks in the jet/outflow and/or to ambipolar diffusion at the wind base. 
The SO abundance in the disk is 3 - 4 orders of magnitude larger than in evolved protoplanetary disks. This may be due to an SO 
enhancement in the accretion shock at the envelope-disk interface or in spiral shocks if the disk is partly gravitationally unstable. 
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1. Introduction 

first steps of the formation of a low-mass star are regulated 
by the simultaneous effects of the mass accretion onto the star 
and the ejection of matter from the stellar-disk system. As a re¬ 
sult, protostars (the so-called Class 0 objects) are characterised 
by the occurrence of fast bipolar jets flowing perpendicular to the 
plane of accretion disks. In practice, although the precise launch 
region (star, inner disk edge at ~ 0.1 AU, outer disk at ~ 0.1- 
10 AU) remains unknown (e.g., iFerreira et al.ll200^ . jets are 
thought to remove excess angular momentum from the star-disk 
system, thus allowing disk accretion onto the central object. Un¬ 
fortunately, the observations of the jet-disk pristine systems in 
deeply embedded protostars are very difflcult to perform given 
the small scales involved as well as the occurrence of numer¬ 
ous other kinematical components involved in the star formation 


recipe (cavities of swept-up material, infalling envelope, static 
ambient cloud). 


The HH 212 region in Orion (at 450 pc) can be considered as 
an ideal laboratory to investigate the interplay of infall, outflow 
and rotation in the earliest evolutionary phases of the star form¬ 
ing process. HH 212 is low-mass Class 0 source driving a sym¬ 
metric and bipolar jet extensively obser ved in typical molecula r 
tracers such as H 2 , SiO, and CO (e.g., IZinnecker et al.l[T998l) . 
High-spatial resolution observations (down to ^ 073 - 0V4) per- 
formed with the SubMillimeter Array (SMA) (iLee et al.ll2006i 
l2007al l2008t). the IRAM Plateau de Bure (PdB) interferome- 


K 


and the At- 


ter (ICodella etaP 120071 l^abrit et al.l l2(M 1^ _ 

acama Large Millim eter Array (ALMA; see lUee et al.l 1201 4t 
ICodella et al.ll2014^ reveal the inner +1" - 2" = 450 - 900 AU 
collimated jet (width ^ 100 AU) close to the protostar. HH 212 
is also associated with a flattened rotating envelope in the equa- 
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Fig. 1. Channel maps of SO 98 - 87 (top) and SO 2 82,6 - (bottom) at systemic, low (LV), intermediate (IV), and high (HV) velocities (panels 
from left to right). The green, blue, and red contours trace the emission at systemic velocity and over the labeled blue- and red-shifted velocity 
intervals (IV - Vsysl). The lowest contour of the emission at systemic velocity is shown in all panels. The first contour is at 5 <t for SO and 3 <t for 
SO 2 with steps of 5cr (VsyO and 20cr (LV, IV, HV) for SO and 3 cr for SO 2 (Iff ~ 3-4 mJy/beam/0.43 kms^*). The tilted black cross is centred at 
the source position and indicates the jet and the disk PA (PAjet = 22°, PAjisi^ = 112°). The ellipse in the bottom-left corner shows the beam HPBW 
of the line emission maps (0'.'65 x 0'.'47, PA = 49°). The triangles indicate the emission peaks along the blue and red lobes. 


tor perpendicular to the jet axis observed firstly with the NRAO 
Very L arge Array (VLA) in NH 3 emission by IWiseman et al.l 
(1200 Ih on 6000 AU scales. More recent SMA and ALMA ob¬ 
servations in the CO isotopologues and HCO^ on ~ 2000, 800 
AU scales indicates that the flattened envelope is not only ro¬ 
tating but also infalling onto the central source, and can there¬ 
fore be identified as a pseudo-disk accord ing to magnetized core 
collapse models dLee et alJ l2006i l2014l) . In addition, a com- 
pact (< 120 AU), o ptically thick dust pea k is observed by 
ICodella et al.l (l2007h : Lee et al.l (l2006l 1201 4l) and attributed to 
an edge-on disk rather than the inner envelope. This seems to be 
confirmed by HCO^ and C'^O emission showing signatures of 
a compact disk of radius ~ 90 AU keplerian rotating around a 
source of ^ 0.2 - 0.3 M© (iLee et alJl2014l:ICodella et^l2014ah . 
Keplerian rotating disks had previously been observed toward s 
only other three Class 0 o bjects: IRAS 4A2 (IChoi et al.ll2010ll . 
LI527 dTobin et al.1 l2012t ISakai et al.1 Imi) . and VLA1623A 
dMurillo et alj|2013l) but HH 212 can be considered as the only 
object clearly revealing both a disk and a fast collimated jet, 
calling for further observations aimed to characterise its inner 
regions. 

Recent observations show that SO can be used to image high- 
velo cify protostellar jets, similarly to a s tandard tracer such a s 
SiO dLee et al.1 1^1 Ot lTafalla et ani2010l: l^della et al.l|2014bl) . 

However, only two protostellar jets have been so far clearly 
mapped in SO lines: HH 211 and NGC1333-IRAS2A, and es¬ 
timates of the SO abundance have been derived only for HH 
211. On the other hand, sulfur bearing species have been long 
searched in protoplanetary disks associated with evolved young 
stellar objects (i.e. Class I-II) but o nly CS has been rou tinely ob¬ 
served while SO is hard to detect: iDutrev et al.l d201 ih reported 


no detectio n of SO (and H 7 S) fr om three prototypical protoplan¬ 
etary disks, iFuente et al.l d201f)l) reported a detec t ion in the disk 
of the T Tauri star AB Aur, and iGuilloteau et ^ d2013h showed 
that SO is exceptionally observed in disks with only one definite 
detection on a sample of 42 T Tauri and Herbig Ae stars towards 
04302-1-2247. The statistics regarding disks around Class 0 pro¬ 
tostars is_even_poorer^he only case is represented by L1527, 
where ISakai et ^ (l2014l) have shown that SO as observed by 
ALMA originates from the outer disk near the centrifugal bar¬ 
rier, and its emission is argued to be enhanced by an accretion 
shock. Observatio ns of HH 212 in t he SO 9% - 87 line obtained 
with the SMA by iLee et al] (l2007al) suggest that SO may have 
two components, a low-velocity one originating in the inner ro¬ 
tating envelope/pesudo-disk and an high-velocity one from the 
jet. However, a detailed analysis of these two components was 
prevented by the lack of angular resolution and sensitivity. 

In this paper we exploit the unprecedented combination of 
high-spatial resolution and high-sensitivity of ALMA to image 
and characterise both the molecular jet and the disk around the 
HH 212 protostar through SO and SO 2 lines. In order to evaluate 
the reliability of SO and SO 2 lines as tracers of shock chem¬ 
istry and/or high-density gas in the disk, the SO and SO 2 spatio- 
kinematical properties are compared with a well-know shock 
chemistry tracer (SiO), with a universal outflow tracer not sensi¬ 
tive to density or chemistry (CO), and with a disk tracer (C^^O). 
The ALMA observations of SO and SO 2 emission presented in 
Se ct. [Hare at higher angular resolution than previous SMA ones 
bv iLeeet al.1 (l2007ah (HPBW = 0765 x 0747, i.e. around a factor 
4 better than at SMA, HPBW = 1716 x 0784) and ~ 100 times 
more sensitive (cr ~ 3 - 4 mJy/beam/0.43 kms“^ with ALMA 
and cr ~ 450 mJy/beam/kms“* with SMA). This allows us to 
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Table 1. Properties of the observed transitions. 


Line 

Vo 

(MHz) 

Enp 

(K) 

(D 2 ) 

logioAij 

(s-') 

SO 2 82,6 - 7ij 

334673.352 

43 

4.9 

-3.9 

SO lOii - lOio 

336553.346 

143 

0.3 

-5.2 

C '^0 3 - 2 

337061.130 

32 

0.04 

-5.6 

CO 3 - 2 

345795.990 

33 

0.04 

-5.6 

SO 93 - 87 

346528.481 

79 

21.5 

-3.3 

SiO 8 - 7 

347330.631 

75 

77.0 

-2.7 


“from the JPL molecular database (Pickett et ah! 1998^ 


disentangle for the first time the origin of the different velocity 
components in SO and SO 2 and to probe the different structures 
in the compact circumstellar region, i.e. the cavity walls, the out¬ 
flow, the molecular jet, the envelope, and the disk (see Sect. O. 
The simultaneous observation of the SiO 8-7 and CO/C'^O 3-2 
lines presented in ICodella et al.l (12014al) allows determining the 
molecules abundances and the physical and dynamical proper¬ 
ties of the jet and the disk (see Sect.|4]l. Finally, our conclusions 
are summarized in Sect. |5] 

2. Observations and data reduction 

HH 212 was observed in Band 7 in the extended configuration 
of the ALMA Early Science Cycle 0 operations on December 
1 2012 using 24 antennas of 12-m. The shortest baseline was 
about 20 m and the longest one 360 m, hence the maximum un- 
filtred scale is of 3" at 850/tm. The properties of the observed 
SO Og-Sv, SO lOii - lOio, and SO 2 82,6 - \,i transitions are 
summarized in Table [T] (frequency vq in MHz, upper level en¬ 
ergy Eup in K, S ij in D^, coefficient for radiative decay Aij in 
s '). We also report the properties of the SiO 8-7, CO 3-2, 
and C"'0 3 - 2 transitions observ ed during the same run and 
presented by ICodella et al.l (l2014al) . which are also analysed to 
compare with the SO and SO 2 emission. The obtained datacubes 
have a spectral resolution of 488 kHz (0.42 - 0.43 km s“ *), a typ¬ 
ical beam FWHM of 0765 x 0747 at PA~ 49°, and an rms noise 
of ~ 3 - 4 mJy/beam in the 0.43 km s“' channel. The calibration 
was carried out following standard procedures and using quasars 
10538^40, J0607-085, as well as Callisto and Ganymede. Spec¬ 
tral line imaging and data analysis were performed using the 
CASA0 and the GILDAl^ packages. Offsets are given with re¬ 
spect to the MMl p rotostar position as det ermined from the dust 
continuum peak by ICodella et al.l (l2014ah . i.e. q'(J 2000) = 05'' 
43“ 51L41, 5(J2000) = -01° 02' 53717. These values are in ex¬ 
cellent agreement with those determined through previous SMA 
and PdBI observations. Velocities are given with respect to the 
systemic velo city, Vsvs, which is es timated following the same 
method as in ICodella et al.l (l2014ah . They show that the emis¬ 
sion in the C'^O 3-2 and C^^'S 7-6 lines is most extended and 
most symmetric in the 0.43 km s“' wide velocity bin centered at 
-1-1.13 ± 0.22 kms^' and -1-1.42 + 0.22 kms^', respectively, and 
assume that Vsys is the average of these values. For SO 98-87 and 
SO 2 82,6 - 7i j the largest extension is observed at -1-1.27 + 0.21 
km s“* and + 1.53+0.22 km s“'. The average of the central veloc¬ 
ity of C'^O, C^^S, SO, and SO 2 gives K v. ^ +1.3 ± 0.2 kms-' , 
in agreement with the value adopted by ICodella et al.l (12014al) . 
Note that despite the coarse velocity binning, all the lines in¬ 
dicate a systemic velocity significantly smaller than what deter- 

' http://casa.nrao.edu 
^ http://www.iram.fr/IRAMFR/GILDAS 


mined bv lWiseman et alJ(l 200 lh f+1 . 6+0.1 kms *) and the value 
adopted bv iLee et all (l200dl2007d IMl (+1.7 + 0.1 kms-'). 
As argued by ICodella et al.l (l2014ah this discrepancy can be due 
to the fact that the C'^O, C^^'S, SO, and SO 2 emission seen by 
ALMA probes the gas motion on much smalle r spatial scales 

£ _ 1350 AU) than the ammonia observed by IWiseman et al.l 

2001 1) (~ 6000 AU). Moreover, SO and SO 2 line peaks towards 
the MMl protostar are redshifted by —(-1 kms-' with respect 
to the velocity where they show the maximum spatial extent (see 
Sect. 13.4b . This may due to optical depth effets, and suggests that 
adopting the velocity where the ambient gas show the maximum 
extent may be more accurate than using the line peak velocity to 
define Vsys. Note, however, that the results presented in the paper 
are not dependent on the 0.3 km s ' velocity difference between 
the different Vsys estimates. 

3. Results 

In order to constrain the origin of the detected SO Os - 87 and 
SO 2 82,6 - 7 i ,7 lines we define four velocity intervals over which 
the emission show different morphologies and kinematic proper¬ 
tied 

- systemic velocity: 0 < |V - Vsysl < 0.6 kms“'; 

- low-velocity (LV hereafter): 0.6 < |y - Vsysl < 5.2 kms-'; 

- intermediate-velocity (IV hereafter); 5.2 < |y - ysysl < 10.2 
kms-'; 

- high-velocity (HV hereafter): 10.2 < |y - y^yj < 15.2 
kms-'. 

Figure [T] shows channel maps of the the SO 93-87 and SO 2 
82,6 “ j emission in the four velocity intervals defined above. 
In the following sections we discuss the origin and the properties 
of the different SO and SO 2 velocity components. 

3.1. The cavity walls 

At systemic velocity SO 93-87 and SO 2 82,6 - lij peak at the 
source position (see left panel of Figure[TJ. This may be due to a 
density enhancement in the circumstel lar region, as sug gested by 
the modeling of HCO^ 4-3 emission (iLee et alj2014ll and/or to 
an enhancement of the SO and SO 2 abundances due to dust grain 
mantles sublimation and release of sulphur-bearing species in 
gas-phase. The SO 93-87 emission also extends in a wide-angle 
biconical structure around the direction of the H 2 /SiO jet which 
extends up to 2" - 3" distance from source and up to +0.6 
km s-' with respect to systemic velocity. Also the SO 2 emission 
shows a similar morphology even though only southern to the 
source. Given the high critical density of the observed SO and 
SO 2 transitions (n^r ~ 2 x 10® cm-^ for gas temperatures of 
50 - 100 K) and t he simi larity with the C^''S emission detected 
by ICodella et al.l (l2014al) . this emission is believed to origi¬ 
nate in the compressed, swept-up gas in the outflow cavity walls. 


3.2. The outflow 

At low velocities the SO 93-87 and SO 2 82,6 - fij emission is 
bipolar and collimated along the jet direction. Blue- and red- 
shifted LV emission larg ely overlap in the two lobe s similarly 
to SiO 5 - 4 and 8 - 7 (ICodella et alJl2(M l2014at iLee et all 
l2007ah . The emission has a transverse FWHM of ~ 0756, which 

^ the defined velocity intervals are the same as in ICabrit et all (l2007l) 
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Fig. 2. From left to right: Position-velocity diagram of CO 3-2, SiO 8-7, SO 98 - 87 , and SO 2 82,6 -7i j lines obtained along the jet (PAjet = 22°, 
top) and the disk (PAdisk = 112°, bottom). Horizonthal and vertical dashed lines mark = -1-1.3 kms^' and the continuum peak MMl. The 
red arrows on the SO 2 82,6 - 'll.! PV indicates emission from complex organic molecules which is blended with the SO 2 emission at the source 
position. For CO 3-2, SiO 8-7, SO 98 - 87 , SO 2 82,6 - lij the first contour is at 5cr (cr = 10 mjy for SiO, 4 mjy for SO, and 3 mjy for SO 2 ) 
with steps of lOcr (CO, SiO, and SO), and 5o- (SO 2 ). 


implies an intrinsic width of ~ 073 = 135 AU after correction 
for the ALMA HPBW (~ 0747 in the trans verse direction), i.e 
larger than the SiO jet width measured by ICabrit et"^ (l2007l) 
(~ 90 AU for all velocity components). This indicates that the 
LV component, with large blue/red overlap, is probing a slower 
wider-angle outflow surrounding the narrow HV jet (see the fol¬ 
lowing section). 

3.3. The fast and collimated molecular jet 

Figure [T] shows that as we move to intermediate and high 
velocities SO Og - 87 and SO 2 82,6 - 7i 7 blue- and red-shifted 
emission only slightly overlaps and become more and more 
collimated. The FWHM is ~ 075 which translate in an intrinsic 
width of ~ 072, i.e. ~ 90 AU, after correction for the beam 
HPBW across the jet (~ 0747 ). Assuming an inclination angle 
i ~ 4° dClaussen et al.irT99^ the de-projected gas velocity is 
V ~ 70 - 210 kms k The inferred velocity and width are in 
perfec t agreement with those estimated from S iO 5-4 emission 
(see, ICodella et alJl 200 ^ ICabrit et alJl2007h and indicate that 
for radial velocities larger than 5 kms * SO and SO 2 trace the 
molecular jet. 

To further investigate the effectiveness of SO and SO 2 as 
tracers of the molecular jet, position-velocity diagrams (PVs) of 


SO and SO 2 emission along and perpendicular to the jet direc¬ 
tion (PAjet = 22 °, PAdisk = 112 °) are extracted and compared 
with the PVs of CO 3-2 and SiO 8-7. Figure |2] shows that the 
considered tracers have different spatio-kinematical properties. 
While SiO traces only jet emission with no contamination from 
ambient gas, SO and SO 2 are dominated by cloud and outflow 
emission at systemic and low velocities. Moreover, the LV and 
IV SO emission shows accelerating “arms”, from systemic ve¬ 
locity at the source position to —h/ - 8 kms^' at ~ +1" - 175 
distance from source, which are not seen in SiO. These arms 
suggest that some SO/S 02 -rich gas from the circumstellar en¬ 
velope or disk may be continuously accelerated to high veloc¬ 
ities away from the source, either by interaction with the fast 
jet or by ma gneto-centrifugal forces ope r ating on envelope/dis k 
scales (e.g., Ciardi & HennebeU3l2oirt iPanoglou et aljboT^ . 
Also CO 3-2 traces multiple components depending on veloc¬ 
ity (the envelope, the outflow, the jet). Moreover, the CO PVs 
show a strong absorption feature at sligthly red-shifted velocity 
(~ -I-1 km s'). Despite the different behaviours of the considered 
tracers, Figure [3] shows that for velocities larger than —h/ - 8 
km s ' the PV diagram of SO along the jet is very similar to 
that of CO 3-2 and SiO 8-7 suggesting that at high velocities all 
tracers probe the same jet component. This in turn allows us to 
compare the emission from SO and CO at high velocities and to 


Article number, page 4 ofM 









































L. Podio et al.: The jet and the disk of the HH 212 low-mass protostar imaged hy ALMA: SO and SO 2 emission 


Table 2. Properties of the observed emission lines on-source (MMl) and at the position of the knots along the jet. 



rms 

T/ a,* 

r max,b 

V 

^ max,r 

FWZI' 

J,* T/ C,* 

r peak 

y d 

^ mb,peak 

jT^bdV^ 

HVjet/disk'^ 

j T^bdV (HVjet/disk)® 


(K) 

(kms^') 

(kms^') 

(kms“ 

') (kms^') 

(K) 

(Kkrns-*) 

(kms^') 

(Kkms-') 

CO 3 - 2 

B2 

0.3 

-16.8 

13 

24.1 

- 

33.9 + 0.3 

370.0+ 1.1 

[-8,-15] 

82.6 + 0.5 

B3 

0.3 

-16.8 

1.1 

24.5 

- 

73.9 + 0.3 

732.5 + 0.9 

[-8,-15] 

46.7 + 0.4 

MMl 

0.3 

-12.9 

7.3 

20.2 

- 

62.9 + 0.3 

424.3+ 1.0 

- 

- 

R2 

0.5 

-9.5 

16.8 

26.2 

- 

71.1+0.5 

567.5+ 1.6 

[+10,+15] 

43.3 + 0.7 

SiO 8 - 7 

B2 

0.5 

-16.8 

8.6 

25.4 

- 2.6 

27.6 + 0.5 

350.2+ 1.6 

[-8,-15] 

87.0 + 0.7 

B3 

0.4 

-15.9 

11.2 

27.1 

-4.7 

25.2 + 0.4 

294.2+ 1.3 

[-8,-15] 

30.5 + 0.6 

MMl 

0.3 

-19.4 

15.9 

35.3 

3.4 

11.2 + 0.3 

163.2+ 1.1 

- 

- 

R2 

0.4 

-7.3 

14.2 

21.5 

9.0 

31.8 + 0.4 

252.6+ 1.1 

[+10,+15] 

32.2 + 0.5 

SO 98 - 87 

B2 

0.1 

-16.3 

6.9 

23.2 

- 8.2 

4.4+ 0.1 

49.4 + 0.4 

[-8,-15] 

18.5 + 0.2 

B3 

0.1 

-15.1 

9.9 

24.9 

-6.4 

13.1+0.1 

140.2 + 0.5 

[-8,-15] 

14.6 + 0.2 

MMl 

0.1 

-12.9 

11.6 

24.5 

1.1 

26.0 + 0.1 

238.1 + 0.4 

- 

- 

R2 

0.1 

-10.3 

12.9 

23.2 

7.7 

11.6 + 0.1 

118.8 + 0.4 

[+10,+15] 

2.2 + 0.2 

SO 2 82,6 - 7 1,7 

B2 

0.1 

-14.1 

1.8 

15.9 

- 10.2 

0.4+ 0.1 

3.7 + 0.2 

[- 8 ,- 10 ] 

0.9+ 0.1 

B3 

0.1 

-11.5 

10.9 

22.4 

-6.5 

1 . 2 + 0.1 

11.6 + 0.3 

[- 8 ,- 10 ] 

0 . 8 + 0.1 

MMl 

0.1 

- 6.0 

11.3 

17.3 

1.0 

4.4+ 0.1 

- 

- 

- 

R2 

0.1 

-7.2 

10.0 

17.2 

- 1.2 

0 . 6 + 0.1 

5.9 + 0.2 

[+ 8 ,+ 10 ] 

0.9+ 0.1 

SO lOii - lOlo 

MMl 

0.1 

-5.2 

4.7 

9.9 

- 

0.7+ 0.1 

4.2 + 0.3 

[-1.9,-3.5] 

1.4+ 0.1 









[+1.9,+3.5] 

0 . 8 + 0.1 

C‘'03 -2 

MMl 

0.1 

-7.7 

10.8 

18.5 

-0.4 

7.7+ 0.1 

30.7 + 0.3 

[-1.9,-3.5] 

5.0+ 0.1 









[+1.9,+3.5] 

3.2+ 0.1 


* velocities are in terms of V - Vsys and the error on the velocity values is half of the spectral element, i.e. ±0.2 km s“* 

“ maximum blue- and red-shifted velocity dehned as the velocity where the intensity becomes lower than lx rms noise 
full width at zero intensity FWZI = [V'max.r - V'max.bl 

velocity at the emission peak. No value is reported for lines affected by self-absorption, i.e. CO 3-2 and SO lOn - lOm 
main-beam temperature at the emission peak. For lines affected by self-absorption r,„b,peak could be a lower limit 
'■ line intensity integrated over the whole velocity profile (i.e. over the FWZI) 

1 high velocity range where jet and disk emission are detected 

^ line intensity integrated over HVjet/disk- For SO lOn - lOio and C'^O 3 - 2, f (HVdisk) is measured for both the blue- and the red-shifted 

line wings 


estimate the abundances of SO and SO 2 in the jet by comparing 
their column densities with that of CO 3-2 (see Sect. 14.2b . 

Figure [T] shows that the SO 98 - 87 and SO 2 82,6 - HV 
emission has three emission peaks, or knots, marked by black 
triangles on the SO and SO 2 channel maps, two along the blue 
lobe (B3 and B2) and one along the red lobe (R2). The knots 
are located at 0'.'68, 1'.'55, and 1"43 respectively, i.e. at a de- 
projected distance from the MMl protostar of around 300, 700, 
and 650 AU. The B2 and R2 SO/SO 2 knots are roughly coinci- 
dent with the B2, R 2 knots observed in SiO 5-4 emission by 
ICodella et al.l (l2007h and th e SN, SS knots in SiO 8-7 and CO 
3-2 emission observed bv iLee et akl (l2007al) . Id spectra of SO 
98-87 and SO 2 82 , 6 - 71,7 are extracted at the position of the knots 
using the same synthetized beam and the obtained intensity pro- 
hles (Tnib in K versus {V - Vsys) in kms^') are shown in Figure 
|4] From the spectra we estimate the emission line properties in 
the knots, i.e.: rms noise in K, maximum blue- and red-shifted 
velocity (Vmax.b and Vmax.r in kms^'), full width at zero inten¬ 
sity {FWZI in kms^'), velocity at the emission peak (Vpeak in 
kms^'), main-beam temperature peak (Tmb.peak in K), and inte¬ 
grated intensity (f T^\,dV in K km s“') (see Table|2]l. 


From the velocity at the SO emission peak reported in Table 
|2]we estimate the knot de-projected velocity, V - VpeakIsin{i), 
and derive the knots dynamical timescale, Tdyn ~ 15 years for 
B3 and ~ 30 years for the symmetrically located blue- and red- 
shifted knots B2 and R2. The position, distance, velocity and 
dynamical timescale of the detected SO knots are summarized 
in Table ID The inferred dynamical timescales are very low sug¬ 
gesting that the observed SO emission probes very recent ejec¬ 
tion events. Note, however, that these estimates are affected by a 
large uncertainty due to the assumed inclination angle, the choice 
of the tracer (Vpeak varies up to a factor 3 depending on the tracer, 
either SO, SO 2 , or SiO), and the broad line prohles (up to ~ 27 
kms^'). The line broadening may be caused by the jet geom¬ 
etry and/or by internal jet shocks. In the first case, if the jet is 
conical and propagates with constant velocity V at an angle i 
to the plane of the sky, the observed maximum blue- and red- 
shifted velocities (see Table |2]i imply a jet half-opening angle 
0m-dx ~ lOfl This in turn implies a knot deprojected velocity 


the jet half-opening angle is tan = (R - 1)/(R + 1) x tan(i), where 
R is the algebrical ratio of maximum to minimum radial velocities and 
i the inclination with respect to the plane of the sky (see ICodella et all 
I2007h . 
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Fig. 3. Left panels: The PV of SO 98 - 87 emission along the jet PA (in 
red) is overplotted on SiO 8-7 and CO 3-2 PVs (in black, top and 
bottom panels, respectively). Right panels: Ratio of the SO 98 - 87 PV 
to the SiO 8-7 and CO 3-2 PV (top and bottom panels, respectively). 
The colour scale is logarithmic. 


V = VmaxI sin(i - 0 tnax) which is about a factor 2 smaller than 
the deprojected Vpeak, implying a factor 2 larger timescales than 
those reported in Table 01 Note, however, that the observed line 
profiles are strongly asymmetric and show broad line wings, un- 
like the sinthetic profile computed bv lKwan & Tademarul(ll98^ 
for a conical jet of constant speed seen close to edge-on (see 
their Fig. 1, dashed curve for i = 82.8° and 20niax = 20°). Al¬ 
ternatively, the observed line broadening may be due to an “in¬ 
ternal jet shock” where fast jet material catches up and shocks 
slower previous ejecta. In this scenario, over-pressured shocked 
material is squeezed sideways out of the jet beam, producing 
an expanding bowshock structure in its wake, which cause the 
observed line broade ning. The bowshock scenario is favored by 
ICodella et al.l (|2007|) and furt h er sup ported by the SiO maps and 
transverse PVs of ILee et al.l (|200^ . In that case, Vpeak traces 
the bow wings while the knot velocity along the jet axis is 

V - Vmixlsin{i), i.e. around a factor 2 larger than the depro¬ 
jected Vpeak- This would imply a factor two smaller timescales 
than those reported in Table 0] We conclude that the estimated 
knot timescales are affected by an uncertainty of a factor 2-3 
depending on the origin of the line broadening and the selected 
knot tracer. 

In Figure0]the obtained SO 98 - 87 and SO 2 82,6 - 7i j inten¬ 
sity profiles are also compared with the profiles of a typic al jet 
tracer, i.e. SiO 8-7 dCodella et aI.ll200A iLee et al.ll2007ah . and 
a jet/outflow (at high/low velocities respectively) tracer, i.e. CO 
3-2, extracted at the same positions and with very similar syn- 
thetized beam (differences between the beams are < O'.'l). The 
comparison shows that even though the SiO and SO knots are 
roughly co-spatial SO and SO 2 peak at higher velocities than SiO 
in the blue lobe and at lower velocities in the red one. This strat¬ 
ification in velocity is a typical signature of shock-excited emis¬ 
sion as both the abundances of the observed molecular species 
and the gas physical conditions (density and temperature) un¬ 
dergo strong gradients in the post-shock cooling region, giving 
rise to chemical segregation. The different profiles and molec- 



Fig. 4. Intensity profiles of SiO 8-7 (blue), SO 98 - 87 (red), and SO 2 
82,6 - 7 1,7 (orange) over CO 3-2 (grey) at the knots (B2, B3, R2) and 
source (MMl) positions. The name and distance of the knots and the 
normalization factor to the CO 3-2 peak intensity are indicated in the 
top of the panels. The orange vertical arrows at the MMl position indi¬ 
cate emission lines from complex organic molecules which are blended 
to the SO 2 line. 


ular stratification shown by the symmetrically located knots B2 
and R2 suggest different shock conditions in the blue and red 
lobes within 2" from source. This asymmetry is in contrast with 
the sy mmetry observed on larger scales in H 2 (IZinnecker et al.l 
1 19981) . and suggests that asymmetries may wash out on large 
timescales. A detailed comparison of the observed intensity pro¬ 
files with shock models predictions is out of the scope of this 
paper and will be presented in Gusdorf et al. (in preparation). 

Finally, it should be noted that, even though SO, SO 2 , SiO, 
and CO show different behaviours at low velocity (i.e. different 
spatio-kinematical distribution in the PVs and different line pro¬ 
files and peak velocities), the maximum blue and red velocities 
are remarkably similar. This indicates that, as also suggested by 
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Fig. 5. Blue- and red- shifted emission of SO lOn - lOio at low velocities (|P - Vsysl < 1.9 kms^', left panel) and in the higher velocity range 
where the emission is consistent with a rotating disk (HVdisk : L9 <\V - Psysl < 3.5 kms^*, middle panel) is compared with C*’0 3-2 emission 
in the HVdisk range (right pane/. l^odella et ^2014ah . The tilted black cross indicates the jet (PA = 22°) and the disk (PA = 112°) direction. The 
ellipse in the bottom-left comer shows the beam HPBW of the line emission maps (0'.'65 x 0'.'47, PA = 49°). First contour at 5cr with steps of 5cr 
for C*’0 3-2 and at 3cr with steps of Icr for SO lOu - lOjo. 


the comparison of the line PVs in Figured all the considered 
lines probe the same jet component at high velocities. 

3 . 4 . The flattened rotating envelope 

As suggested by iLee et al.1 (12007 ah while the high-velocity SO 
Og-Sv emission is clearly associated with the molecular jet, the 
low velocity one could arise in the inner rotating envelope and/or 
in a pesudo-disk. The bottom panels of Fig. |2] shows the PVs of 
the observed emission lines obtained cutting the datacube per¬ 
pendicular to the jet PA, i.e. along the disk direction (PAdisk = 
112°). With the exception of SiO, which shows poor emission 
along the disk, the other lines show a velocity gradient perpen¬ 
dicular to the jet PA of a few km s ' over a scale < 1" (i.e. 
<450 AU), which suggests rotation in a flattened envelope and/or 
a disk with the same rotation sense as observed in HCO^ by 
iLee et all (|2014 and in C'^O 3 - 2 bv ICodella et'aD(l2014ah . 

The SO Og -87 and SO 2 82 , 6 -7i j intensity prohles at the po¬ 
sition of the MMl protostar (Fig.|4]i, peak at sligthly redshifted 
velocity, i.e. ~ -i-l kms *, where CO 3-2 shows a deep ab¬ 
sorption feature. This could be a signature of infalling gas. We 
also note that the SO 2 82,6 - 7i 7 line profile shows two spectral 
“bumps” at blue- and red-shifted velocities, which are not de¬ 
tected in the other tracers. The PV diagram in Fig.|2]shows that 
these spectral features are seen only in a spatially unresolved 
region towards the MMl position, and not along the outflow. 
This suggests that SO 2 is blended with two unidentified lines 
emitted in the inner few 100 AU of the protostellar envelope. 
These lines are likely produced by Complex Organic Molecules 
(COMs) evaporated from the icy mantles of dust grains in the 
warm circumstellar region (Tdust >100 K), as observed toward s 
many other low-luminosity protostars (e.g., iMaurv et al.ll2014l) . 
The identification and analysis of these lines is out of the scope 
of this paper. 

3.5. SO IO 11 -IO 10 : the compact rotating disk ? 

Channel maps and position-velocity diagrams along and perpen¬ 
dicular to the jet direction are also obtained for the fainter SO 
lOii - lOio line as shown in Figures|5]and|6l At difference with 
SO 9 g -87 and SO 2 82,6 - 7ij, the emission in the SO lOn-lOio 
line is compact (~ 90 AU scales) and narrow in velocity extend¬ 
ing only up to +5 kms^' with respect to systemic. As shown 
in Fig. |5] the red- and blue-shifted emission peaks are displaced 
roughly along the jet direction for velocities <1.9 kms * (see 


left panel). However, in the velocity range where C'^O 3-2 
probes a compact (~ 90 AU) disk r otating around a 0.3 ± 0.1 M© 
star (1.9 < |y - Vsysl < 3.5 kms '. ICodella et ^l2014ah the SO 
lOii-lOio emission peaks are consistently displaced along the 
disk axis (see middle and right panels). This suggests that the 
SO emission at those velocities may also originate in the com¬ 
pact disk. 

The PV diagram perpendicular to the jet axis in Fig. sup¬ 
ports this scenario as it shows a blue- and a red- lobe which 
are consistent with emission from a rotating structure. However, 
when cutting the datacube along the jet PA the blue- and red- 
shifted lobes are consistent with jet emission. This indicates that 
SO lOii-lOio traces partly rotating gas in a flattened envelope 
and/or in a disk, and partly the jet. 

In Figure|2]the SO lOn-lOio line intensity profile extracted 
at the source position is compared with C'^O 3-2 analysed by 
ICodella et al.l (1201 4ah . In the 1.9 < |y - ysyj < 3.5 velocity 
range the SO 10n-10io and C'^O 3-2 line profiles are very 
similar suggesting that at those velocities they are tracing the 
same disk component. However, at low velocities the SO lOn- 
lOio line intensity profile shows a red-shifted absorption feature 
at ~ 0.5 kms '. This may be due to the SO lOn-lOio higher 
excitation energy which makes it less affected by the emission 
from the circumstellar gas and more sensitive to absorption by 
the infalling gas. As the observed SO lOn-lOio is compact, the 
detected absorption feature could probe the infall of the outer 
disk regions, on much smaller scales than the infalling gas traced 
by CO 3 - 2. 


4. Discussion 


4.1. Searching for jet rotation 

Magneto hydro-dynamical (MHD) models predict that the jet is 
launched by magneto-centrifugal forces which lift the gas from 
the disk along the open magnetic held line, accelerating and col¬ 
limating it. According to these models the jet rotates around its 
axis carrying away the excess angular momentum from the disk, 
thus allowing accretion onto the central star. Many recent studies 
have been devoted to search for systematic velocity asymmetries 
across the jet axis as a signature of jet rotation, first at op tical 
wavelengths (e.g., iBacciotti et al.ll200^ICoffev et al. 2004) and 
more recently in the su b- mm and mm range (e.g., Codgl Ia_et al.l 


I20071 II .ee et al .mOOTafl l 2 ^ l2009l llTunhardt et al .Il 200 ^ 

Ten tative rotation mea s urements have been o btain ed for HH 
212 by iDavis et ^ (l2000h . ICodella et al.l (l2007h . and iLee et al.l 
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SO 10,1-10,0 - PA=22 SO 10,,-10,0 - PA=112 



Fig. 6. From left to right: Position-velocity diagram of SO lO,, - 10,0 
obtained along the jet (PAje, = 22 °, left) and the disk (PAjisi^ = 112 °, 
right). Horizonthal and vertical dashed lines mark Vsys = -t-1.3 kms“* 
and the continuum peak MMl. The first contour is at 6 cr (Itr = 2 mjy) 
with steps of O.Str. On the PV along the disk the keplerian curves for 
M, = 0.3 ±0.1 Mq (solid and dashed black curves) and the r“* curve 
for infalling gas with angular momentum conservation (red curve) are 
overplotted. 



-10 -50 5 10 

V — V (km/s) 


Fig. 7. Intensity profiles of SO 10,, - 10,0 (grey) and C'^O 3-2 (red) 
extracted at the source position. The horizonthal bars in the bottom of 
the panel indicate the high velocity range (HVdis k = 1-9 - 3.5 kms~‘ ) 
where the lines are consistent with a disk origin dCodella et alj|20143) . 
The normalization factor to the SO 10,, - lOio intensity at peak (solid 
line) and in the blue-shifted HVdisk range (dotted line) are labeled on the 
top of the panel. 


(l2008h . iDavis et ^ (l2000l) reported a rotation speed of ~ 1.5 
km s * at a distance of ~ 230 AU from the jet axis in the south¬ 
ern H 2 knot SKI located at 2300 AU distance from the source. 
However, the velocity asymmetry in the northern lobe would in¬ 
dicate rotation in the opposite sense. Moreover, at such large dis¬ 
tances from the driving source and from the jet axis the velocity 
pattern is likely dominated by effects other than rotation, i.e. by 
interaction with the surrounding cloud in t he bow-shock wings , 
and/or jet wiggling and precession (e.g., ICorreia et al.ll 200 ^ . 
More recent high angular resolution measurements in the inner 
B2 and R2 knots (~ 1" - 2 " from source) gives contradictory 
results: ICodella et akl (|2007|) do not observe transverse velocity 
shifts above 1 kms“* in the SiO 5 - 4 emission with 0778 x 0734 
HPBW resolution, while higher angular resolution SMA obser¬ 
vations of SiO 8-7 (HPBW = 0736 x 0733) indicate a velocity 
shift at the tip of the two bow-shocks (iLee et al.ll2008h . 

To test the validity of these tentative rotation measurements 
we search for velocity asymmetries across the collimated jet by 
extracting PV diagrams of SO Og-Sy and SO 2 82,6 -7, 7 emission 



0.5 0 - 0.5 0.5 0 - 0.5 

offset from jet axis (") offset from jet axis (") 


Fig. 8. From left to right: Position-velocity diagrams of SO 98 - 87 , 
and SO 2 82,6 - 7,7 lines obtained perpendicular to the jet axis (PA = 
112°) at the position of knots B2, B3, and R2. The distance of the knots 
from the driving MMl protostar is indicated. The first contour is at 5cr 
(cr = 4 mJy for SO, and 1.5 mJy for SO 2 ) with steps of lOcr (SO), and 
5o- (SO 2 ). 


perpendicular to the jet axis at the positions of knots B2, B3, and 
R2 (see Fig.[ 8 ]l. PV diagrams of SiO are n ot shown a s these are at 
lower resolution with respect to those bv iLee et alJ (l2008ll . The 
figure shows no evidence of a velocity gradient across the jet 
axis, consistent with the high-velocity SO/SO 2 jet being practi¬ 
cally unresolved transversally. Higher angular resolution (~ 071) 
which will be available with the full ALMA array is needed to 
properly sample the velocity pattern across the jet width and to 
provide a reliable test of rotation in the HH 212 jet. 

In the PV across the B3 knot there is a hint of rotation for the 
gas close to systemic velocity, which could indicate rotatio n of 
the outflow cavity as seen in (see lCodella et al.ll2014ah . 


4.2. SO and SO 2 abundance in the jet 


In order to constrain the physical conditions of the gas and the 
SO abundance in the inner jet knots (B3, B2, and R2) the ob¬ 
served SO line intensities (corrected for beam dilution) and ra¬ 
tios (independent of beam dilution) are compared with the pre¬ 
dictions of the statistical equilibrium, one-dimensional radia¬ 
tive transfer code RADEX a dopting plane parallel slab geom¬ 
etry dvan der Tak et al.ll2007h . Figure^ shows the SO line tem¬ 
perature ratio, Tr(SO lOn - 10 io)/Tr(SO 98 - 87 ), versus the 
SO 98 - 87 intrinsic line temperature, Tr(SO 98 - 87 ), pre¬ 
dicted by RADEX for a reasonable range of kinetic temperatures 
(Tk = 50 - 500 K ) and Hy densities (n/y, from 10^ to 10^ cm“^) 
(e.g., ICabrit et"^ 2007h and for SO column densities (Nso/AV) 
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Table 3. Beam-averaged column densities and abundances of SO, SO 2 , and CO in the jet knots B2, B3, and R2. The column densities are 
estimated from the SO 98 - 87 , SO 2 82,6 — lij, and CO 3-2 line intensities integrated over the indicated high velocity ranges assuming optically 
thin thermalized emission at = 50 - 500 K. The abundances of SO and SO 2 are derived as ^specie = A^specie/^co x Alco/A^m and assuming 
A:co = NcoINn, = lO-^ 


Jet 

knot 

HVje, 

(kms-') 

Nso 

( 10 '"' cm-®) 

Nco 

(lO'® cm-®) 

Xso 

( 10 -®) 

HVjet 

(kms-') 

Nso2 

(lO'® cm-®) 

Nco 

(lO'® cm-®) 

Xso2 

( 10 -®) 

B2 

[-8,-15] 

3.0-7.1 

>0.4- 2 . 1 * 

<7.5-3.3* 

[- 8 ,- 10 ] 

0.1 - 1.3 

> 1.8 - 10 . 1 * 

< 4.7 - 12.4* 

B3 

[-8,-15] 

2.3 - 5.6 

0 . 2 - 1.2 

10.5-4.6 

[- 8 ,- 10 ] 

O.l - 1.1 

p 

OA 

1 

p 

4.7-12.4 

R2 

[-h10,-h15] 

0.4 - 0.8 

0 . 2 - l.l 

1 

p 

bo 

[-h 8 ,-h 10 ] 

0.1 - 1.3 

- 

- 


* As CO 3-2 is optically thick in knot B2, the estimated CO column density, Nco, is a lower limit, and the SO and SO 2 abundances, X^q and 3fso,, 
are upper limits. 
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Fig. 9. SO line temperature ratio Tr(SO lOn - 10 io)/Tr(SO 98 - 87) 
versus intrinsic line temperature Tr(SO 98 - 87) predicted by RADEX 
slab models at Tr = 50, 500 K (black and red points/curves). Each curve 
corresponds to the labelled column density NsolXV increasing from 
10*^ to 10'^ cm“^ (kms^*)^' by factors of 10 (left to right). Each dot 
along the curves corresponds to the labeled logfu^j) (from 10 ^ to lO’ 
cm“’). Optically thin models are marked by a cross, optically thick one 
by a star. The black points indicate the upper limit of the SO line ratio 
and the SO 98 - 87 intensity corrected for beam dilution in the B2, B3 
and R2 knots over the HVjet velocity range. 


increasing from 10'^ to 10^^ (kms~^)“'. The model predic¬ 

tions are compared with the observed line ratios and line intensi¬ 
ties cotTected for beam dilution. As discussed in Sect. l3.3l for ve¬ 
locities higher than ~ 8 kms“^ with respect to systemic the SO 
emission is co-spatial with that of CO 3-2 and SiO 8-7 indicat¬ 
ing that SO originate in the collimated molecular jet with no con¬ 
tamination from the surrounding outflow and cloud. Therefore, 
the SO ratios and intensity are measured on the jet high-velocity 
range defined as: HVjet = [- 8 , -15] km s * for knots B3 and B2 
and HVjet = [- 1 -10, H-15] km s ' for knot R2. The different veloc¬ 
ity range defined for R2 is to avoid the strong absorption feature 
around -1-9 km s ' when comparin g wit h the CO 3-2 emis sion. 
As discussed by iLee et all (12007 ah and ICabrit et'aP (1201 2h this 
feature is present at all positions and may be caused by an ex¬ 
tended foreground component fully resolved-out by interferom¬ 


eters. The SO lOn - lOio line is not detected in the jet knots, 
hence we report the upper limit of the considered SO line ratio. 
The SO line ratios do not depend on the assumed beam Ailing 
factor, as long as both lines trace the same volume of gas. The 
observed SO 98 - 87 line intensity, T^b, is corrected for beam 
dilution, as the jet of ~ 072 width is broadened by a factor ~ 2 
by beam convolution across the jet. The errorbars on the SO line 
intensity in Fig.j^accounts for the range of intensity values over 
HVjet, and for possible beam dilut ion along the jet ax is (a factor 
~ 3 if the knot is circular, see also lCabrit et alJl2007h . 

Figure 0 shows that the avalaible SO observations do not 
allow constraining the density and temperature of the emitting 
gas. However, for the considered temperatures and H 2 densities 
the observations are in agreement with the model predictions for 
column densities NsoIXV ~ 10'^ - 10'"' cm“^ (kms^')^'. For 
these values of the column density the emission in the consid¬ 
ered SO lines is optically thin. 

If the density is sufficiently high, i.e. for n/fj larger than a 
few 10 ® cm^^, the emission is also in local termodynamic equi¬ 
librium (LTE). In this case, an estimate of the beam-averaged 
SO column density in the inner jet knots can be derived from 
the SO 98 - 87 line intensity integrated over the high velocities 
range HVjet assuming optically thin and thermalized em ission at 
Tk ~ Tex = 50 - 500 K. Following ICabrit et alJ (1201 2l) also the 
beam-averaged CO column density is derived from the CO 3-2 
intensity integrated on the same velocity range as for SO assum¬ 
ing optically thin, LTE emission. The latter assumption is not 
valid for knot B2 where the ratio between the CO 3-2 and SiO 
8-7 emission is ~ 1 over the HVjet velocity range. This indicates 
that both lines are optically thick. Therefore, for knot B2 the de¬ 
rived CO column density is a lower limit. Then, the abundance of 
SO is derived as Xso = Xco x Nso /Nco, where the abundance of 
CO with respect to H 2 is assumed to be Xco = A^co/Mij = 10“"'. 
The main uncertainty on the estimated column densities and 
abundances is due to the assumption on the gas temperature, 
therefore we report a range of values corresponding to the tem¬ 
perature values Tjc ~ Tex = 50 - 500 K. Moreover, if the density 
in the jet is lower than a few 10® cm“^, the SO Og - 87 line is 
sub-thermally excited implying higher SO column densities. In 
particular, assuming nH 2 = 3 x 10® cm ®, the SO column densi¬ 
ties and abundances estimated using RADEX are from a factor 
2 (for Tk = 500 K) to a factor 4 (for Tk = 50 K) larger than in 
the LTE-optically thin case. 

The SO 2 82,6- 7 1,7 emission is much fainter and covers lower 
velocities than SO, therefore it is not possible to compare its PV 
with that of CO 3-2 (see Eig. |2]). However, based on the inten¬ 
sity profiles shown in Fig.|4]we tentatively assume that SO 2 and 
CO 3-2 emission probe the same jet component over the veloc¬ 
ity range HVjet = [-8, -10] kms“' in knots B3 and B2 and that 
both are optically thin and thermalized. Based on these assump- 
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tions the column density and abundance of SO 2 are inferred by 
applying the same method as for SO. The main source of error 
on the estimated SO 2 column densities and abundances is due to 
the assumed gas temperature, therefore we report a range of val¬ 
ues for Tk ~ Tex = 50 - 500 K. As for SO, if the density in the 
jet is lower than a few 10® cm“^, the SO 2 line is sub-thermally 
excited. In that case, using RAD EX with = 3 x 10® cm ® 
and Tex = 50 - 500 K we obtain a factor 2 to 4 higher column 
densities and abundances. 

The SO, SO 2 , and CO line intensities integrated over the high 
velocities range HVjet in the inner jet knots (B3, B2, and R2) are 
reported in the last column of Tab.|2l while the column densities 
and abundances estimated assuming LTE-optically thin emission 
are summarized in Tab. [3] The SO abundance in the high veloc¬ 
ity jet of HH 212 (Xso ~ 10“^ - 2 x 10“®) is com parable to 
that e stimated in other Class 0 jets, namely HH 211 (iLee et akl 
l2010ll . the extremely high velocity (E HV) wings of the L 1448- 
mm and IRAS 04166-1-2706 outflows (iTafalla et al.ll2010h , and 
the B1 and B2 bowshocks at the tip of the LI 157 outflow cavity 
dBachiller & Perez Gutierrez! 19971) . On the other hand, the ratio 
SO 2 /SO appears variable from jet to jet: we find S 02 /SO~ 0.5-4 
in HH 212, while it is around 1 in LI 157, and < 0.1 - 0.3 in the 
L1448-mm and IRAS 04166-H270 EHV gas. 

Assuming a sulfur elemental abundance (S/H)o ~ 1.38x10“® 
(lAsplund et alJ l2005h . these values indicate that from 1% up 
to 40% of the elemental S is in the form of SO and SO 2 in 
the HH 212 jet. This represents an enhancement by a factor 
10-100 with respect to SO and SO 2 abundances in prestel- 
lar cores and protostellar envelopes dBachiller & Perez Gutierre^ 
ll997HTafallaetalJl2010h . 

SO and SO 2 enhancements up to ^ 10“^ - 10 ® are 

predicted by models of magnetized molecular C-shocks, as 
a result of neutral-neutral re actions in the shock-heated gas 
dPineau des Eorets et al.l Il9 93[) and sputtering of H 2 S-rich icy 
grain mantles dElower & Pineau des Foretsi[l994t) . The forma¬ 
tion of SO 2 by oxydation of SO is predicted to occur late in the 
post-shock, hence the ratio of SO 2 /SO would be an indicator o f 
shock age (see e.g. Fig 7 of iFlower & Pineau des Foretslll994l) . 
In the case of HH 212, the ratio SO 2 /SO > 0.5 that we observe 
seems inconsistent with the very short dynamical timescale of 
15-30 years of the knots. Such a high ratio at early times would 
require that SO 2 is mainly released from grain mantles rather 
than formed in the gas phase out of SO as assumed in published 
shock models. Observations of other sulfur-bearing species, such 
as CS, OCS, H 2 S, H 2 CS, are crucial to constrain the sulfur chem¬ 
istry and obtain an accurat e estimate of the su lphur released in 
gas-phase in shocks (e.g., iPodio et^l2014t) . Alternatively to 
shocks, high SO and SO 2 abundances up to 10“® are predicted at 
the turbulent in terface between the outflow and the surrounding 
cloud material dViti et al.ll2002h . SO and SO 2 is also strongly en¬ 
hanced in an MHD disk wind due to the thermal sublimation of 
icy mantles near the source and the grad ual heating by ambipo - 
lar diffusion during MHD acceleration (iPanoglou et al.llT012h . 
Detailed comparison with such models will be the subject of a 
future paper (Tabone et al., in prep). 

4.3. Jet physical and dynamical properties: H 2 density and 
mass loss rate 

From the beam-averaged CO column density estimated in the 
HVjet range, Nco (see Tab. O, one can infer the H 2 column 
and volume density in the knots along the jet, Ah, and hh,, and 
the jet mass loss rate, Mjet. The determination of these quanti¬ 
ties relies on the assumption that CO 3-2 emission at high ve¬ 


locity is optically thin and traces all the emitting mass in the 
beam. As CO 3-2 in knot B2 is optically thick the Ah,, wh,, 
and Mjet values inferred for this knot are lower limits. First, the 
beam-averaged H 2 column density is derived assuming Xco = 
Aco/A^h, = 10“^*. Then, the H 2 density is derived assuming that 
the jet is a cylinder uniformely filled by the observed gas. Hence, 
mh, = (7/PBWtr/2/?jet) X A^H,/ 2 /?jet, where 2Rjet is the jet width 
(~ 90 AU) and the factor (HPBWal2Rjet) accounts for the beam 
dilution of the ~ 90 AU jet emission in the transverse beam di¬ 
rection. Finally, the mass loss rate is estimated assuming that the 
mass in the jet flows at constant density and speed along the jet 
axis over the beam length. The latter assumption is not fulfilled 
if the gas in the kno t s is highl y compressed by shocks, therefore 
following iLee et af] (l2007albh we correct the mass loss rate for a 
compression factor of ~ 3. Based on the above assumptions, Mjet 
is calculated as Mjet = 1/3 x mn, x (Nco/Xco) x HPBWa x Vjet, 
where Vjet is the de-projected jet velocity in the considered HVjet 
range (Vjet ~ 165 kms“') and HPBWt, is the beam size across 
the jet width {PfPBWi^ - 0747). Note that the correction for 
compression is uncertain, as it depends on the unknown shock 
parameters (magnetic held and shock speed), therefore the Mjet 
values are affected by a factor ~ 3 uncertainty. Moreover, the 
estimated mass loss rate accounts only for the mass transported 
by the molecular jet component at high velocity (|y - Vsysl > 8 
kms“'). The total jet mass loss rate may be up to a factor of a 
few higher than the estimated value as the low velocity compo¬ 
nent can considerably contribute to the total mass loss rate if it 
traces a slow w ide angle wind rath er than entrained ambient ma¬ 
terial (see, e.g.. lMaurri et al.ir2014t) . If the jet is partially ionized 
the atomic gas may further contribute to the mass loss rate. How¬ 
ever, Spitzer and Herschel observations of molecular jets driven 
by Class 0 sources indicate that the mass flux transported by the 
atomic component is about 10 times smaller than that carried 
by th e molecular gas (iDionatos et al.l 1200^ 120101: iNisini et aH 
l2015h . Hence, the atomic component of HH 212, if any, should 
not signiflcatly contribute to the total jet mass loss rate. 

The inferred values of Ah,, uh,, and Mjet for the B2, B3, and 
R2 knots are summarized in Tab.|?| These quantities are affected 
by the uncertainty on the Aco estimates, hence a range of values 
is given depending on the temperature assumed to derive Aco 
(Tex = 50,500 K). As explained above, the Mjet values inferred 
for the molecular high velocity jet component are a lower limit to 
the total jet mass loss rate. Therefore, the mass flux transported 
by the jet (Mjet ^ 0.2 - 2 x 10“® Mq yr“') is at least 3% to 33% 
of the envelope infall rate (Mint ~ 6 x 10 ® M^ vr“^ iLee et all 
l2006h . indicating an high jet efficiency {MytIMmf ^ 0.03 - 0.3), 
in agreement with the young age of the source. Note that the 
mass accretion rate from the disk onto the central source may 
be higher than the envelope infall rate. However, there are no 
estimates of the disk mass accretion rate for HH 212 as Class 0 
sources like this one are too deeply embedded in their parental 
envelope to enable the use of well-calibrated accretion tracers 
such as, e.g., the Bry line. 

4.4. Disk properties: mass and SO abundance 

An analysis similar to that performed in Sect. I4.2l is applied to 
constrain the gas physical conditions and the SO abundance in 
the disk. The observed SO lOn - lOio line intensity and the 
SO lOii - lOio/SO 98 - 87 ratio in the blue and red disk lobes 
(HVdisk ; |V - Vsysl = 1.9 - 3.5 km s“^ Ic^lla et al.llM4al) 
are compared with the predictions of RADEX computed for 
Tk = 50 - 500 K, H 2 densities larger than 3x10® cm“® and SO 


Article number, page 10 ofM 














































L. Podio et al.: The jet and the disk of the HH 212 low-mass protostar imaged hy ALMA: SO and SO 2 emission 

Table 4. Properties of the knots along the HH 212 jet: position offset (RA and dec in "), on-sky and deprojected distance from the MMl protostar 
(dtan in " and d in AU), deprojected peak velocity (V in km s“‘), dynamical timescale (Tdyn in yr), high velocity range on which the jet properties 
are estimated (HVjet in kms^'), beam-averaged H 9 column density (Nn^ in cm“^), H 2 volume density («h 2 in cm^'^j, and mass loss rate (Mjet in 
Moyr-‘). 
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RA 
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^tan 
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(AU) 

V 
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^dyn 

(yr) 
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B3 
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1 
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(-/I 
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-1.43 
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30 
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0.2- 1.1 X lO^' 

0.3- 1.6 X 10® 

0.2- l.Ox 10-® 


* As CO 3-2 is optically thick in knot B2, the estimated CO column density, Nco, and hence H 2 column and volume density, and and 
mass loss rate, Mjet, are lower limits. 
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Table 5. Beam-averaged column densities and abundances of SO and 
C*’0 in the disk. The column densities are derived by integrating the 
line intens ity of SO lOn - lOip and C*’0 3-2 over the velocity range 
defined bv lOidella et al.l ( l2014iJl (HVhm.: |V-Kysl = 1.9 - 3.5 kms-‘). 
We assume optically thin emission at LTE with Tex = 50 - 500 K. The 
abundance of SO is derived as X, 

A^co/ A^Ho and assuming Nr nn/Xco = 1/1792 and Xco 
lO '* dWilson & Roodlll^ . 


specie A^specie/A^Cl’O ^ ^ 


A^co/A^h, = 



Disk 

HVdi.sk 

Nso 

Nquo 

Aso 

- 

lobe 

(kms) 

(lO'® cm-^) 

(10^® cm-2) 

(10-^) 

- 

blue 

[-1.9,-3.5] 

6.0-4.6 

0.2- 1.3 

1.4-0.2 


red 

[-h1.9,-h3.5] 

3.4-2.6 

0.2 - 0.9 

1.2-0.2 


Fig. 10. SO line temperature ratio Tr(SO 10ii-10io)/Tr(SO Og-Sv) ver¬ 
sus intrinsic line temperature Tr(SO IOh - lOio) predicted by RADEX 
slab models at Tr = 50, 500 K (black and red points/curves). Each 
curve corresponds to the labelled column density Nso/XV increasing 
from 10*^ to 2 X 10*^ cm“^ (kms^*)^' going from left to right. Each 
dot along the curves corresponds to the labeled log(n//j) (from 10 ^'^ to 
> 10* cm“*). Models for which the emission in the SO lOn - lOio is 
optically thin are marked by a cross, optically thick ones by a star. The 
black points indicate the SO line ratio and the SO lOn - lOio line inten¬ 
sity (or their lower/upper limits) corrected for beam dilution in the B2, 
B3 and R2 knots and in the disk blue and red lobes. 


column densities Aso/AV = 10'^,5x 10'^, 10^®, 2 x lO'® 
(kms^')^'. The observed SO lOn - lOio line intensity, Tnib, is 
corrected for beam dilution, as the emission from the compact 
disk lobes (/?disk ~ 90 AU, i.e. 0'/2) is strongly diluted in the 
0'/65 X 0'.'47 beam (dilution factor ~ 8 ). The SO % - 87 emission 
at low velocities is heavily dominated by the outflow emission, 
hence the observed SO lOn - lOm/SO Og - 87 ratio is a lower 
limit to the SO line ratio in the disk. Figure [TO] shows that, to 
reproduce the bright SO lOn - lOm emission, column densi¬ 
ties of Aso/AV ~ 1 - 2 X 10 '® cm“^ (kms^')^' are required, 
i.e. two orders of magnitude larger than what estimated in the 
jet (Aso/AV ~ 10 '"^). In the considered models the emission in 
the SO 98 - 87 line is optically thick (see also Fig.|9]), while the 
emission in the SO lOn - lOio line is still optically thin. More¬ 


over, volum e densities are e xpected to be larger than 10^ cm“^ 
in the disk (iLee et al.ll2014ll . therefore we can assume that the 
line emission is thermalized. 

As the emission in the SO lOn-lOio is optically thin and 
thermalized, we estimate the SO column density in the disk by 
integrating the SO lOn-lOio line intensity on the velocity inter¬ 
val where the emission is consistent with a disk origin, ITVdisk, 
and assuming Tk ~ Trx = 50 - 500 K. Give n its low abundance 
(Acno = Arn/1792. IWilson & Roodlll994h and critical density 
the emission in the C'^O 3-2 line is also assumed to be opti¬ 
cally thin and thermalized at Tk ~ Tex = 50 - 500 K and the 
C'^O column density is inferred from the C'^O 3-2 line inten¬ 
sity integrated on ITVdisk as for SO. The SO lOn-lOio and C'^O 
3-2 line intensity integrated on HVdisk are reported in the last 
column of Table |2] while the estimated beam-averaged SO and 
C'^O column densities are summarized in Table|5] 

From the beam-averaged C'^O column density, Nqiiq, we 
estimate the H 2 column and volume density in the disk and 
its mass. The beam-averaged H 2 column density is derived as¬ 
suming Xq\iq - Aco/ 1792 and Aco = Aco/A^Hj = 10“"^. 
Then, as the disk is seen almost edge-on, the H 2 density is 
8 X Ah 2 //^disk. where /?disk ~ 90 AU and the factor 8 accounts 
for the beam dilution. Finally, the disk mass is obtained by inte¬ 
grating the beam-averaged C*^0 column density over the beam 
area. Abeam, and summing over the blue and red disk lobes: 

Afdisk ~ (1.4 X TOHj X(Aci7o/Acl7o)xAbeam)b|ue+red’ where the 

factor 1.4 accounts for the mass in the form of Helium. 

Depending on the adopted temperature value, the estimated 
beam-averaged H 2 column density and H 2 volume density in the 
disk are Ah 2 ~ 0.4-2 x 10^^ cm“^ and «h 2 ~ 0.2 - 1 x 10^ cm“^, 
where the lower and higher values correspond to Tgx = 50, 500 
K, respectively. The disk mass turns out to be ~ 0.002 - 0.013 
M q. This value is c lose to the disk mass of 0.014 Mq estimated 
bv lUee et al.l (l2014t) from HCO^ 4-3. Hence, if C'^O is tracing 
the same disk layer as HCO^, it is indeed optically thin and all 
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in the gas phase in this region. Note that more mass could be 
hidden in the dark disk mid plane where molecules are frozen 
onto dust grains, so our estimate is only a lower limit to the total 
disk mass, i.e. Mdisk > 0.002 - 0.013 Mq. 


Finally, we estimate the SO gas-phase abundance in the disk 
by assuming that SO lOn-lOio and C'^O 3-2 originate from the 
same disk layer in the 1.9 - 3.5 kms“^ velocity range and that 
C'^O is all in gas-phase in the considered disk-emitting region. 
Based on the above assumptions Xso = x Nso/Ncno ~ 

10“^ - 10^^. This estimate would remain correct even if C'^O 
is depleted onto grains, as long as SO and C^^O are depleted by 
the same amount (and both optically thin). The estimated SO col¬ 
umn density and abundance are much higher than what observed 
in evolved protoplanetary disks (Nsct < lO'^ cm“^, Xso < 10“^', 
iFuente et al.ll2010l:lDutrev et alj|2()l Ih and may be due to an SO 
enhancement in the a ccretion shock a t the envelop e -disk inter- 
face as suggested by iLee et all (l2014l) : ISakai et all (I20l4l . Al¬ 
ternatively, an high SO abundance (up to 10 is predicted by 
self-gravitating disk models due to the shocks occuring in the 
disk spirals caused b y gravitational instabilities (lllee et alJl201 lb 
iDouglas et al.ll20f3h! Since the central star has a m ass of only 
0.2 - 0.3 Mo Lee et akllml ICodella etaDl2014ah and our es¬ 
timates of the disk mass do not include the mass reservoir frozen 
onto grains, such a scenario is not excluded. Flowever, the pre¬ 
sented ALMA observations do not allow us to distinguish be¬ 
tween these scenarios. The angular resolution (~ 0'.'5) is too low 
to fully resolve the disk, therefore the SO 10n-10io line profile is 
averaged over a large portion of the disk impeding to detect SO 
enhancements as a function of velocity. Moreover, an edge-on 
view is not favorable to search for SO enhancements in a grav¬ 
itationally unstable disk as each line of sight would cross many 
spiral arms. Such a scenario could however be tested for less in¬ 
clined Class 0 disks with higher resolution ALMA observations. 


5. Conclusions 


HH 212 is the prototype of the Class 0 jet and it has been the sub¬ 
ject of a number of high angular resolution stud ies at mm wave¬ 
lengths perfo rmed with the SMA and the PdBI (Lee et aL|[2006l 
l2007al 120081:ICodella et aLll2Q07l: ICabrit et alJl2007ll2012ll . The 


high sensitivity of the ALMA Cycle 0 observations allow us for 
the first time to image not only the envelope, the outflow, and 
the jet but also the in ner compact disk (see also lLee et al.ll2014l : 
ICodella et al.ll201^ . In particular the presented analysis of SO 
and SO 2 emission show that these molecules probe different 
components in the inner ~ 1300 AU of the HH 212 protostellar 
system : (i) the circumstellar gas and the X-shaped cavity walls 
at systemic velocity; (ii) the poorly collimated gas (width ~ 135 
AU) likely swept out by the high-velocity molecular jet at low 
velocities; (iii) the collimated (~ 90 AU) and fast (~ 100 - 200 
kms^') molecular jet at intermediate and high velocities; (iv) a 
compact structure with the blue- and the red-shifted peaks dis¬ 
placed roughly on the equatorial plane, which is consistent with 
an origin in the 90 AU dis k rotating a round a ~ 0.2 - 0.3 Mp 
star pre viously detected by iLee et all (l2014ll and ICodella et al.l 
(I2014ah . in the high-excitation SO lOn - lOio line. 

The high-velocity molecular jet detected in SO, SO 2 , CO, 
and SiO lines is dense (~ 10^ - 10^ cm“^) and transport a mass 
loss rate Mjet > 0.2 - 2 x 10“® Mq yr“', indicating a high ejection 
efficienc y (M\alM m S 0.03 - 0.3, for an infall rate of ~ 6 x 10“^ 
Moyr \ iLee et all l2006h . The abundance of SO in the jet is 
Xso ~ 10^^ - 10“®, similar to what estimated in other Class 


0 outflows and jets (e.g.. 

Bachiller & Perez GutierrezI 119971: 

iTafalla et al.ll2010tlLee et alJ 

2 OIOI). while Xsn,/Xso > 0.5. This 


indicates that between 1% and 40% of the elemental S has been 
converted into SO and SO 2 in the shocks occuring along the jet 
or at the outflow-cloud interface d Pineau de s Forets et al.lll993l ; 

I Flower & Pineau des Foretsill994l: IViti et al. 20021) and^ at the 


base of the wind due to ambipolar diffusion ( Panoglou et al.l 
| 2 ^ . 


From C'^O 3-2 emission we estimate a disk mass > 
0.002 - 0.013 Mq, which is con sistent with the estimate derived 
from HCO^ bv iLee et al.l (l2014l) . The SO abundance in the disk 
(Xso ~ 10“^ - 10“^) is 3 - 4 orders of magnitud e larger than 
what estimated in evolv ed protoplanetary disks (iFuente et al.l 
l2010l:[Dutrev et al.l201 Ih . Such an SO enhancement may be pro¬ 
duced in the accretion shock occuring at the interface betw een 
the envelope and the disk as suggested bv iLee et al.l (1201 4l) and 
ISakai etal.1 (12014 . or in spiral shocks if the disk is part ly gravi¬ 
tationally unstable (lllee et al.l201 ItlPouglas et al.l2013h . Higher 
angular resolution and sensitivity observations are required to 
address the exact origin of the detected high-excitation SO emis¬ 
sion. 
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